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Abstract 

Gate operations in a quantum information processor are generally realized by tailoring specific 
periods of free and driven evolution of a quantum system. Unwanted environmental noise, which 
may in principle be distinct during these two periods, acts to decohere the system and increase the 
gate error rate. While there has been significant progress characterizing noise processes during free 
evolution, the corresponding driven-evolution case is more challenging as the noise being probed 
is also extant during the characterization protocol. Here we demonstrate the noise spectroscopy 
(0.1 — 200 MHz) of a superconducting flux qubit during driven evolution by using a robust spin¬ 
locking pulse sequence to measure relaxation ( T \ p ) in the rotating frame. In the case of flux noise, 
we resolve spectral features due to coherent fluctuators, and further identify a signature of the 
1 MHz defect in a time-domain spin-echo experiment. The driven-evolution noise spectroscopy 
complements free-evolution methods, enabling the means to characterize and distinguish various 
noise processes relevant for universal quantum control. 
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Introduction 


Within the gate-based model of quantum computation, a quantum algorithm may be 
efficiently implemented on a quantum information processor by decomposing it into a finite 
and discrete set of coherent gate operations [1]. Their realization, although in detail depen¬ 
dent on the particular physical-qubit modality, may be broadly divided into two categories, 
free-evolution and driven-evolution, indicating the period(s) during which the desired gate 
operates on the system. During these two types of evolution, unwanted environmental noise 
acts to decohere the system and increase the gate error rate. While process tomography 
may provide the complete map of the action of a gate operation, it does not characterize the 
noise sources which cause fidelity degradation and which, in principle, may differ under free- 
and driven-evolution conditions. In contrast, the noise sensitivity of the quantum device 
itself, in conjunction with tailored pulse sequences, may be used to identify spectral features 
of the noise and thereby elucidate the underlying noise mechanisms [2, 3]. 

Dynamical decoupling protocols, e.g., spin echo and its multi-pulse extensions, have been 
applied to mitigate dephasing during free evolution in numerous qubit modalities, including 
atomic ensembles [4, 5] and single atoms [6], spin ensembles [7], diamond nitrogen-vacancy 
centres [8, 9], semiconductor quantum dots [10, 11], and superconducting qubits [12], The 
spectral filtering properties [13] of these and related sequences have proven a useful tool for 
performing free-evolution noise spectroscopy with superconducting qubits over a wide range 
of frequencies spanning millihertz [14, 15] to over 20 MHz [12], Several groups have also 
made progress characterizing [12, 16, 17] and mitigating [18] decoherence during driven evo¬ 
lution. This case is more difficult in practice, however, because driven-evolution decoherence 
may be highly correlated with the operational errors that occur while implementing a par¬ 
ticular protocol. Such correlations pose challenges for fault-tolerant quantum information 
processing, and it motivates the importance of characterizing the noise processes which are 
manifest while the qubit is under external drive. 

Noise spectroscopy during driven evolution has been studied most extensively via Rabi 
nutation of the qubit between its ground and excited states, achieved by continuously driving 
the qubit with an oscillating field. Noise power at the Rabi frequency is inferred from the 
Rabi-oscillation-decay envelope, and the noise spectrum is reconstructed over the achievable 
range of Rabi frequencies by varying the field amplitude [12, 16, 17] (and F. Yoshihara et 
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al., manuscript in preparation). There are several drawbacks to this approach, however, 
including the Rabi experiment’s associated transverse decay, its sensitivity to low-frequency 
fluctuations of the Rabi frequency, the interpretation of the resultant non-exponential decay 
law, and the practical consideration of sampling sufficiently many points to resolve and fit 
accurately a decaying sinusoid oscillating at the Rabi frequency. This type of Rabi-based 
noise spectroscopy can, however, be generalized by preparing the qubit in an arbitrary 
initial state. One unique case, originating within the NMR community, is the so-called 
Tip experiment [16, 19-22], which measures the driven-evolution analog to 7\ relaxation. 
In Tip , the driving held is applied along an axis collinear with the qubit’s state-vector in 
the rotating frame, resulting in a condition called “spin locking.” The main advantages of 
Ti p -based noise spectroscopy are that the qubit decays longitudinally from the spin-locked 
state with a straightforward exponential decay law (no oscillations); the decay is dominated 
by noise at the spin-locking frequency and, like its T\ analog, is relatively insensitive to 
broadband low-frequency noise; and it can therefore have higher accuracy over a wider 
frequency range than the standard Rabi-based approach. 

In this work, we develop a modified T\ p pulse sequence to perform noise spectroscopy 
during the driven evolution of a superconducting flux qubit. The modified pulse sequence 
addresses several non-idealities due primarily to 1/f-type noise that may pollute the spin¬ 
locking dynamics when using the conventional T lp sequence. Using this improved sequence, 
we measure both the flux and tunnel-coupling noise spectra during driven evolution over 
the frequency range 0.1 — 200 MHz. The observed flux-noise power spectral density (PSD) 
generally agrees with the 1/f dependence found independently using free-evolution noise 
spectroscopy [12]. Notably, between the experiments in Ref. [12] and those presented here, 
we thermal-cycled the device from 12 mK to 4K and back again. With the higher accuracy 
afforded by the T lp measurement, we resolved two additional “bump”-like features in the 
spectrum, possibly due to a set of two-level systems (TLSs), e.g., electron spins, that were 
apparently activated by the thermal cycling. We further demonstrate that the underlying 
noise mechanisms associated with these particular spectral features are also active during 
free evolution by observing their signature in a time-domain echo experiment. 
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Results 


Device. Our test device is a superconducting flux qubit (Fig. la), an aluminum loop 
interrupted by four A1-A10 Z -A1 Josephson junctions. When an external magnetic flux $ « 
$ 0/2 threads the loop, the qubit’s potential energy assumes a double-well profile. The 
wells are associated with clockwise and counterclockwise circulating currents of magnitude 
I p ~ 180 11 A and energies ±he/2 = ±/ p ( < F — *3?o/2) tunable by the applied flux. These 
diabatic circulating-current states tunnel-couple with a strength A = 5.4 GHz, and the 
resulting hybridized qubit states are manipulated using microwave pulses that couple to 
the qnbit via an on-chip antenna. The qubit is embedded in a DC SQUID, a sensitive 
magnetometer which serves as the qnbit readout. The device and measurement set-up is 
detailed further in Ref. [12]. 

Reference frames and Hamiltonians. In the laboratory frame and within a two-level 
system approximation, the qubit Hamiltonian is: 

i-L = (h/ 2) [(A + 5A(t))a z + (e + 6e(t))i t x + A ri (t) cos(27rz/ rf f + (p)a x ] , (1) 

where a z and & x are Pauli operators. 5A(t) and Se(t ) represent random fluctuations in the 
effective tunnel-coupling and flux terms, respectively, and the oscillating term represents 
the harmonic driving field which is treated semiclassically with amplitude A r f(f), frequency 
zz r f, and phase <fi. Although noise affecting the qubit physically manifests itself in the lab 
frame (Eq. 1), the longitudinal (Tj and Tj p ) and transverse (X 2 and T ^ p ) relaxation times for 
free- and driven evolution may be analogously defined in the qubit-eigenbasis and rotating 
frames, respectively, as described below and presented in Table 1 (also see Supplementary 
Notes for detailed derivation and description). Their measurement may then be related back 
to the physical noise parameters in the laboratory frame. 

Decoherence during free evolution is described by transforming Eq. 1 to the qubit’s energy 
eigenbasis: 

W = (ft/2) [„„ dy + A ri (t) cosd cos(27rz/ rf t + 0)dy 
+ A r f (t) sin 9 cos(27rz/ r ft + (j))a z i 

+ (<5A (t) cos 6 + 8e[t ) sin 9) g z i + (— 5A(t) sin 9 + 5e(t ) cos 9) dy] , (2) 

where 9 = arctan(£/A) is the rotation angle of the quantization axis from the lab frame, 
and z/ q = \Je 2 + A 2 is the qubit’s frequency splitting (Fig. lb). The standard longitudinal 
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(Xi = 1/ri) and transverse (T 2 = 1/^) relaxation times are defined in this frame. In our 
device, T x « 12 /xs with only a weak dependence on e, whereas T 2 ~ 3 /xs at e — 0 and 
decreases away from this point [12]. 

Decoherence during driven evolution (i.e., during a constant-amplitude pulse, A r{ (t) = 
A r f) is described by a second transformation, taking Eq. 2 to a reference frame rotating at 
the drive frequency z/ r f (omitting stochastic terms): 

% = (/i/2) [Az/ a z + z/ R (cos 0ay + sin 0ay)] , (3) 

where A v = z/ q — z/ r f is the frequency detuning between the qubit and the driving field, 
and z/ R = ^A r fCosd is the Rabi frequency under resonant driving (Az/ = 0). Here, the 
rotating wave approximation has been invoked by omitting the fast oscillating terms (see 
Supplementary Notes for a generalization). 

Analogy between free and driven evolution. For resonant driving along the A-axis 
(Az/ = 0, 0 = 0), the Hamiltonian in Eq. 3 describes a pseudospin quantized along X with 
a level splitting equivalent to the Rabi frequency z/ R . The qubit’s Bloch vector precesses 
around the X-axis at a rate z/ R , and its subtended angle depends on the initial state. That 
is, the driven evolution of the qubit state may be interpreted as the “free evolution” of 
a pseudospin around a static field z/ R (Fig. lc). By analogy with the conventional free- 
evolution T\ and T 2 times, we define longitudinal and transverse depolarization times T\ p 
and T 2p for the pseudospin with respect to this new quantization axis (see Table 1). Just as 
T\ depends on the noise at the qubit frequency z/ q that is transverse to the 0-axis (Eq. 2), 
so is Tip determined by the noise at the pseudospin’s level splitting z/ R that is transverse 
to the A"-axis (Eq. 3). For completeness, the measurement and value of the rotating-frame 
decoherence time T 2p (“free-induction” time in the rotating frame) is equivalent to the 
conventional Rabi experiment and associated decay time Tk a bi (Supplementary Fig. S2b). 
More detailed discussion about the analogy between free evolution and driven evolution in 
the rotating frame can be found in Supplementary Notes. 

The connection between the driven-evolution times (T lp , T 2p ), the free-evolution times 
(Tf, T 2 ), and the PSD of the physical noise terms (SA(t), Se(t )) is formally made using the 
generalized Bloch equations (GBE) [23, 24] and the transformations between Eqs. 1-3. The 
exponential relaxation rate during weak and resonant driving is 

Tip = — =-Ti + T^ , (4) 

J- ip ^ 
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where 


= \ Sz{v R ) = \ [cos 2 d SaM + sm 2 d S&J] . (5) 

S\(f) = (2tt) 2 JZ dt|(A(0)A(t) + A(t)A(0)) exp(— i2nft) stands for the symmetrized PSD of 
A-noise, and A G {A, £ ... or x',y',z'...} indicates either the fluctuating part of the system 
Hamiltonian or the axis along which those fluctuations occur. The relaxation rate Tip arises 
from two terms related to noise in the qubit eigenbasis (Eq. 2). The first term (Ti term) is 
associated with transverse noise along the x'-axis at frequency u q ± z/ R (which appears as z/ R 
in the rotating frame). The second term (T^ term) is the longitudinal noise along the Y-axis 
at the drive frequency z/ R , which is directly transferred into the rotating frame. This latter 
term provides a direct method to extract the noise PSD at the locking Rabi frequency by 
measuring relaxation during driven evolution via a T lp experiment. 

Advantages of T\ p over Rabi spectroscopy. The T\ p experiment, when utilized as a 
noise spectrum analyzer to extract the PSD, has inherent advantages over direct Rabi-based 
techniques [12, 17] (and F. Yoshihara et ah, in preparation). One major shortcoming of 
the Rabi method is its vulnerability to low-frequency fluctuation of the effective Rabi 
frequency. First, the Rabi frequency is first-order sensitive to fluctuations of the driving- 
held amplitude Su R (proportional to <T4 rf ), whose RMS fluctuation is found in this case 
proportional to z/ R (<5z/ R ^0.06%z/ R ) [18]. Hence, the temporal z/ R inhomogeneity becomes 
important when the driving is strong (large z/ R ). In addition, low-frequency fluctuation 
of the qubit’s level splitting (5u) modulates the effective Rabi frequency to second order, 
adding to the inhomogeneous broadening [12], The effect dominates the Rabi decay when 
the driving is weak enough (small z/ R , as the second-order contribution (5z/ 2 )/z/ R exceeds 
decay rates from other sources) (Table 1). These two effects are analogous to the scenarios 
of free-evolution dephasing due to low-frequency noise via linear and quadratic coupling, 
respectively [16]. In contrast, the locking dynamics in the Tip experiment largely reduce its 
sensitivity to these fluctuations, and thus substantially improves the accuracy of extracted 
noise PSD, especially at both high and low frequency regimes. 

In addition, T lp has a simple exponential decay law, whereas the Rabi decay law can be 
complicated by its sinusoidal nutation in conjunction with a non-exponential decay envelope 
related to its sensitivity to low-frequency noise (see Table 1), both of which pose a challenge 
to high-precision fitting. A Ti p measurement requires far fewer datapoints than would a 
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Rabi decay envelope on a sinusoid oscillating at the Rabi frequency. These properties make 
the Tip method more accurate and experimentally efficient. 

Conventional T lp sequence (SL-3). The measurement of the rotating-frame relaxation 
time Tip begins by preparing the qubit state to be either parallel or anti-parallel with X 
(collinear with the driving held), corresponding to the pseudospin’s ground and excited 
state, respectively. This results in the so-called “spin-locking” condition, a reference to the 
pseudospin remaining aligned in the absence of noise with the static driving held z/ R in Eq. 3 
(see Fig. le and Supplementary Fig. S2a). The measurement of the relaxation time from 
the spin-locked condition to a depolarized steady-state is the so-called Tj p experiment. 

The original T ip sequence as developed in NMR and applied to superconducting qubits 
was a three-pulse sequence [25], two 7t/2 pulses separated by a 90°-phase-shifted continuous 
driving pulse (labelled here “SL-3”, see Fig. Id). The phase shift is implemented by altering 
the driving-held phase (j) via standard IQ mixing techniques. The dynamics can be visualized 
with the assistance of the Bloch sphere within the rotating frame (Fig. le). For convenience, 
we use {X,Y,Zj to denote rotating frame coordinates, pulse phase, and qubit states. The 
hrst vr/2_ pulse (bar indicates negative or 180-degrees-rotated orientation) brings the qubit 
from its ground state to the equator, parallel with X (step I). The driving held is then applied 
along X and thereby aligned with the qubit state, so that the qubit is effectively locked in 
this orientation and experiences relaxation T lp only (steps II and III). After a finite duration 
r, the remaining polarization along X is projected, by the last 7t/2_ pulse, back along Z 
for readout (step IV). The long-time, steady-state population carries no net X-polarization, 
because the device is operated in the limit z/ R <C k B T/h = 1.3 GHz, where T « 65 mK 
is the effective device temperature obtained by measuring the SQUID’s switching-current 
distribution. The Y-X-Y phase order applied here and the X-Y-X one in [25] are effectively 
the same, since only the relative phase is important. We repeat the same sequence on a 
nominally identically prepared system a few thousand times to improve state estimation. 
Hence, decoherence of this system should be understood in the sense of a temporal ensemble 
average. 

Modified T\ p sequence (interleaved SL-5a and SL-5b). In solid-state systems, the 
recorded T lp decay signal using the original (SL-3) sequence may suffer from several types 
of system imperfection, in particular, signal distortion due to low-frequency noise via other 
channels. One example of a measured T lp trace is shown in Fig. 2b (brown trace). Another 



example can be found in Ref. [25], an early demonstration of spin-locking implemented on 
a quantronium qubit. A prominent feature in these examples is the presence of unwanted 
Rabi-frequency oscillations on the back of an otherwise exponential decay function. We 
excluded the possibility of a constant frequency detuning in our experiment, which could 
effectively separate the driving held and qubit state during the locking period and lead to 
oscillations. Rather, we determined that the oscillations were due to dephasing during the 
intervals between adjacent pulses. In our experiments, we use Gaussian pulses with ~10ns 
width for the 7t/2 pulses and Gaussian rise-and-fall for the long driving pulse. To clearly 
separate pulses, there is a practically unavoidable duration of free evolution between pulses 
(Fig. Id). Depending on the length r of the second pulse, phase diffusion accumulated 
during the first interval will sometimes be doubled (for z/ R r = 1,2,3...), and sometimes 
be refocused (for z/ R r = |, l|,2|...), producing those oscillations. In the presence of this 
noise, the locking protocol is compromised to a differing degree for the temporal ensemble 
elements, making the observed decay a complicated combination of precession, relaxation 
and decoherence processes. To correct for it, we developed a modified five-pulse sequence 
(labelled “SL-5a,” see Fig. 2a), which adds one 7r Y pulse in the middle of each interval in 
SL-3. The added pulses refocus the phase diffusion as is done in a spin-echo sequence [26]. 
Relaxation obtained using SL-5a is shown in Fig. 2b, and it shows a clear improvement (no 
oscillation) over SL-3. 

In several measured traces, we observe a second type of signal distortion—random fluc¬ 
tuations at the few-percent level on the exponential signal, and at a time scale of minutes 
or tens of minutes—e.g., Fig. 2c (blue trace). As we describe in detail below, low-frequency 
fluctuations of the frequency detuning result in the “shivering” of the effective Rabi field 
and, thereby, fluctuations in the observed decay signal. We first describe how we corrected 
the distortion, and then explain its origin. The solution is to use a complementary pulse 
sequence (labelled “SL-5b,” see Fig. 2a) together with SL-5a. The new sequence is nominally 
identical to SL-5a, except that we replace the two 7r/2 _ pulses with two 7t/2 v pulses, such 
that the qubit Bloch vector is anti-parallel with the driving axis for SL-5b. Ideally, both 
sequences should give the same decay function. However, when we apply the sequences in 
an interleaved sampling order (inner loop: alternate between SL-5a and SL-5b; outer loop: 
step r), the anti-symmetry of the signal fluctuation, slow compared with the collection time 
at each value r, is captured by the twin sequences, as shown by the example in Fig. 2c. We 
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then recover the smooth exponential decay by taking the average of the two traces. 

Onr ability to recover a smooth exponential decay is closely related to the ultra-slow 
(~minutes) fluctuation of frequency detuning 8u. First, the occurrence of the signal distor¬ 
tion at large r indicates a deviation of the steady-state population from zero X-polarization. 
Assuming small detuning (Ato first order, the steady-state X-polarization [24] (and 
Supplementary Notes), in the practical limit is R ^k B T/h<^v q , is 

sin r] S x (is q ) 

\ S x {v q ) + S z {uJ ’ 



where r/ = arctan(Au/u R ) is the angle of the driving held with respect to X due to frequency 
detuning (see Fig. 2d). Equation 6 is applicable to both SL-5a and SL-5b, but the difference 
in the last back-projection pulse gives rise to a differential signal response to the non-zero 
X-polarization. Therefore, fluctuation of the frequency detuning, to first order, causes 
fluctuation of the steady-state longitudinal polarization and the anti-symmetric feature of 
the traces in Fig. 2c. The relaxation dynamics in both sequences with finite detuning are 
sketched in the right panel of Fig. 2d. Ignoring the process at small r where the much faster 
T- 2 P process is involved, the normalized decay signal, to first order of rj, can be approximated 
as P a (r) = |(1 — (<T\-) SS )(1 — e^ VlpT ) for SL-5a, and P b (r) = |(1 + (<Jx) ss )(l — e^ VlpT ) for SL- 
5b. The average P(r) = \{P a {r) + P b (r)) = \{1 — e~ VlpT ) becomes independent of detuning, 
and recovers the underlying relaxation signal. 

Noise spectroscopy during driven evolution. We use the modified T\ p pulse se¬ 
quence to implement noise spectroscopy during driven evolution (Fig. 3). To derive the 
A-noise PSD, we perform a driven-evolution Tj p measurement for various u R at £ — 0, 
where the qubit is first-order insensitive to e noise, as well as the standard (free-evolution) 
inversion-recovery Ti-measurement [26]. Both Tj and Tip traces are fit to extract the expo¬ 
nential decay rates Ti and Ti p , which give the A-noise PSD by S a(^ r ) = S z >(v R ) = 2Fi p —Ti. 
For e-noise PSD, we perform the same experiments at e / 0, where the qubit becomes pre¬ 
dominantly sensitive to e noise. We then have S £ {v R ) = [»SA(zqJ — cos 2 6 i S'a(w)]/ s i r i 2 0 = 
[2r lp — Ti — cos 2 6 ) S'a(^' r )]/sin 2 ^, where is previously measured. The locking Rabi 

frequency, z/ R , is determined by measuring Rabi oscillations with the same driving-field 
amplitude. The lower limit (~ 100 kHz) of the spectrum is determined by the Au inhomo¬ 
geneity (Su 2 ), which breaks the locking condition when driving is too weak. On the other 
hand, an undetermined effect (presumably heating due to strong driving) which distorts 
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the SQUID’s switching-current distribution prevents us from probing higher frequency than 
~ 200 MHz. We believe the limitation is unrelated to the SQUID plasma mode, which is 
around ~ 2.3 GHz for this device. 

In this device, the measured A-noise spectrum is generally white. Since the extrapolated 
1// A-noise from low-frequency measurements crosses the white noise level around the lower 
limit, we are unable to determine whether the 1// noise extends to this regime or has a cut¬ 
off [14]. The noise could be due to Johnson noise from charge or critical current fluctuations, 
or clue to thermal photon noise [27] originating from the LC resonator associated with the 
readout SQUID which inductively couples to the qubit (Fig. la). 

The £-noise spectrum exhibits a 1//“ dependence with a = 0.9, a value that gener¬ 
ally agrees with independently measured free-evolution noise spectroscopy using a Ramsey 
(1 mHz — 100 Hz) [14] and CPMG (0.1 — 20 MHz) approach [12]. Although not necessarily 
a universal result, this suggests that driven evolution does not activate e -noise sources for 
this device. In general, the T lp spectroscopy confirms these results with better accuracy and 
extends the frequency range by 5x (presumably limited in this device by a heating effect at 
high driving fields). Note that measuring at higher frequency (0.1 — 1 GHz) is feasible in 
principle by using this technique. 

Spectral and temporal signatures of coherent TLSs. In contrast to previous free- 
evolution measurements, we observe on the 1/f trend two clear “bump-like” features, one 
around 1 MHz and the other around 20 MHz. These features appeared after thermal-cycling 
the device from base temperature to 4.2 K and back again, which likely activated the fluctua- 
tors. The one at 1 MHz has a clear peak above the general 1/f background. This excludes the 
possibility of random telegraph noise (RTN), possibly generated by some microscopic two- 
level fluctuators (TLFs), as the cause of this bump, since RTN would produce a Lorentzian 
spectrum centered at zero frequency [28, 29]. To fit the main features, we assume that the 
e-noise PSD is a combination of the 1/f 0 ' 9 line and two Lorentzians. Both bumps are better 
fit to Lorentzians centered around 1 MHz and 20 MHz than those centered at zero frequency. 
One example is shown by the black line in Fig. 3. 

We find a signature of the 1 MHz bump feature in a free-evolution, time-domain echo 
experiment. Fig. 4 shows the spin-echo phase-decay as measured at several flux-bias points 
where the qubit coherence is mostly limited by £ noise. Under a Gaussian-noise assumption, 
the echo sequence is most sensitive to the noise whose frequency is about the inverse of the 
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total free-evolution time of the sequence. The echo decay functions all exhibit a clear “dip” 
feature around 1 /xs, corresponding to the spectral feature around 1 MHz. To model this 
result, we use the previously assumed PSD function S'(f) (l// 09 plus two Lorentzians) with 
tunable parameters to reproduce the echo-decay signal by the filter function method [12, 
16, 30]. The agreement is good, and bump parameters estimated from the temporal echo 
measurement are within 30% of a direct fit to the driven-evolution spectral data. In general, 
we are able to predict the phase decay from the 7% spectroscopy results, and, in turn, justify 
the performance of the Tj p spectroscopy method on a superconducting qubit. In the case of 
the 20-MHz Lorentzian, simulations indicate that the echo signature corresponding to this 
broad, relatively small feature in the noise spectrum would yield only a small deviation from 
the spin-echo phase-decay function with no discernable ’’dip.” The temporal resolution in 
Fig. 4 is furthermore insufficient to resolve such a deviation. We note that on a subsequent 
cooldown, with the qubit exhibiting a similarly small spectral feature, we did observe a small 
distortion (but no “dip”) in the expected spin-echo phase-decay function. 

It is possible that the observed Lorentzian bumps are generated by two sets of coherent 
TLSs [31], such as electron spins. In a magnetic environment, each spin executes Larmor 
precession at a frequency proportional to its local held. Its transverse component couples to 
the qubit’s flux degree of freedom with a strength depending on geometric relations among 
the device, the held and the spin [32], and adds to a net fluctuating hux whose PSD rehects 
the spectral (Larmor-frequency) distribution of spins, even though the ensemble has already 
lost its coherence. The width of the Lorentzian indicates in part the held inhomogeneity. 
Assuming uncorrelated surface electron spins and an average coupling strength of 2.7 x 
fr° m simulation, both Lorentzian bumps correspond to ~ 10 6 spins. The number 
scales down in the presence of a high degree of spin order [33]. Moreover, if we assume the 
spins feel the same magnetic held as the qubit, about 0.5mT ( B 0 = <3>/Al, where <f> = $ 0/2 
is the hux through the qubit loop and A l ~ 2/irn 2 is the loop size), the corresponding 
electron-spin Larmor frequency is 14 MHz. This is only a rough estimate, as screening due 
to superconducting metal may lead to a large variation of the held at various locations in 
the vicinity of the metal surface, and other locations for these spins, e.g., superconductor 
insulator boundary [34] and substrate, are also possible. Nonetheless, the crude estimate is 
consistent with our observations. 
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Discussion 


The modified T ip method presented here, in conjunction with existing free-evolution 
spectroscopy methods, enables the characterization of the noise sources relevant for both 
driven- and free-evolution quantum control methods. Although many noise sources may be 
manifest similarly during both driving conditions, in principle, the noise in these two cases 
need not be identical. For example, driven evolution may activate certain noise sources 
that are otherwise dormant during free-evolution. In this demonstration with a flux qubit, 
the flux and tunnel-coupling noise spectra reveal more detailed information and cover a 
wider frequency range, a lOx increase in this device, as compared with our previous free- 
evolution (CPMG spectroscopy) and driven evolution (Rabi spectroscopy) methods [12]. We 
observed the spectral features of two sets of coherent TLSs in the environment, possibly due 
to effective electron spins on the metal surface, which are active during driven evolution. We 
could furthermore observe a temporal signature of one of these fluctuators in a free-evolution 
echo experiment, indicating that it is active during both driven- and free evolution. This 
type of noise characterization serves as an important step towards the engineered mitigation 
of decoherence through improved materials, fabrication, control sequences, and qubit design. 
Our modifications to the Tj p pulse sequence target the unwanted effects of low-frequency 
noise, opening this method to the multiple qubit modalities that are subject to such noise. 

Methods 

Description of the qubit. We fabricated our device at NEC, using the standard Dolan 
angle-evaporation deposition process of Al-AlO^-Al on a Si02/Si wafer. 

Our persistent-current qubit [12, 35] consists of a superconducting loop with diameter 
d ~ 2 /i nr, interrupted by four Josephson junctions. Three of the junctions each have the 
Josephson energy Ej — h X 210 GHz, and charging energy Eq = h x 4 GHz; the fourth is 
smaller by a factor a = 0.54. 

The qubit is embedded in a hysteretic DC SQUID, a sensitive magnetometer which 
is used for the qubit readout. It has a critical current / csq = 4.5 /iA; normal resistance 
= 0.25 kf2; mutual qubit-SQUID inductance Mq_ s = 21 pH; on-chip shunt capacitors 
(G s h ~ 10 pF, inductors L sq « 0.1 nH, and bias resistors Ri = lkf2 and Ry = 1 kO. The 
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shunt capacitors bring the plasma frequency down to f p = 2.1GHz, which is much higher 
than the frequencies of the coherent fluctuators presumably related to the frequency- and 
time-domain results. 

Pulse generation and calibration. The r.f. microwave pulses are created by mix¬ 
ing the in-phase and quadrature pulse envelopes, generated by a Tektronix 5014 Arbitrary 
Waveform Generator (AWG), with a continuous wave, provided by an Agilent E8267D PSG 
Vector Signal Generator. We use the internal I/Q mixer of the PSG for the envelope mixing. 
All pulses are further gated in order to eliminate leakage. 

Imperfections in the electronics and coaxial cables will cause pulse distortions, especially 
for short pulses (a few nanoseconds long). To ensure that the pulses we send to the cryostat 
are free from distortions, we first determine the transfer function H ext with a high-speed 
oscilloscope, and then use Hffl to correct for imperfections in the AWG and in the I/Q 
mixers [36]. This set-up allows us to create well-defined Gaussian-shaped microwave pulses 
with pulse widths as short as 2.5 ns. 
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Table 1: Comparison of decohering properties of a flux qubit during free and driven evolution 
(weak and resonant). 
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Figure 1: Device schematic and spin-locking dynamics, (a) Schematic of the flux qubit with 
measurement set-up. (b) Measured qubit frequency versus flux bias, (c) Analogy between free- and 
driven evolution. Red and yellow arrows indicate the quantizing field, respectively. The longitudinal 
axes are highlighted by green arrows, and the transverse planes are highlighted by the blue disks. 
Grey arrows indicate the longitudinal (Xi and T \ p ) and transverse ( T 2 and X 2 p) depolarization, 
(d) Standard three-pulse spin-locking sequence (SL-3). Green and blue indicate X- and V-pulse, 
respectively, (e) Bloch sphere representation of the rotating-frame qubit dynamics under SL-3. The 
purple arrows are the qubit’s state polarization, while the magenta arrows indicate the driving-field 
orientation. The qubit is initially prepared in its ground state (I). The first 7t/2_ pulse rotates the 
qubit by 90° into the equatorial plane (II). The second 90°-phase-shifted continuous driving pulse 
of duration r, is then aligned with the qubit state, effectively locking the qubit along X. During 
the pulse, the qubit undergoes relaxation in this rotating frame towards its steady state (III). The 
final vr/2_ pulse projects the remaining polarization onto Z (=z') for readout (IV). 
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(a) Robust sequences 
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Figure 2: Robust spin-locking sequences, (a) Modified five-pulse sequences (SL-5a and SL- 
5b), robust in the presence of 5e and 5 A low-frequency noise, (b) Comparison of measured signals 
by SL-3 (brown) and SL-5a (blue). Both traces are taken at e = 192 MHz and is R = 3 MHz. Pgw 
is the measured switching probability of the readout SQUID, and is linear with the qubit’s state 
population, (c) Comparison of measured signals by SL-5a (blue), SL-5b (red) and their average 
(purple). Both traces are taken at e = 960 MHz and zx R = 24 MHz. The data are taken in an 
interleaved order between blue and red traces, (d) Cross-section sketch of the qubit dynamics 
under the SL-5a and SL-5b protocols, with the driving field (magenta arrow) tilted by an angle ij 
(r) > 0 in this example) due to slow fluctuations of the frequency detuning. On the left are the 
qubit states right before the locking pulse; on the right, their relaxation paths during the locking 
period are indicated by the dashed lines. The steady-state solution (the contacting point of blue 
and red dash lines) deviates away from {ax} = 0 when r/ ^ 0. 
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Figure 3: Noise spectroscopy during driven evolution, e (red dots) and A (blue dots) noise 
spectroscopy by the Ti P method. Here, previous results on e noise by the CPMG (grey dots) and 
Rabi (green dots) methods measured on the same device are plotted for comparison. From the 
fitting results, the error rate in the worst case is about 30% for e noise, and 70% for A noise. 
Previous experiments and the T\ p experiment were separated by six months and a 4K warm-up. 
The red and blue solid lines are the l/f a power laws, extrapolated from separate low-frequency 
noise measurements [14]. The black line represents the function S'(f) = A/ /°' 9 + S\C(f\ F\, W\) + 
S 2 £(f;F 2 ,W 2 ), where C(f\F,W) = W 2 /((/ - F) 2 + W 2 ) , and A = (2vr0.65 x 10 6 ) 2 (rad-s" 1 ) 2 , 
Si = 7.0 X 10 8 rad-s _1 , F\ = 1.05 MHz, W\ = 0.25 MHz, S 2 = 1.2 x 10 7 rad-s -1 , F 2 = 20MHz, 
W 2 = 25 MHz. 
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Figure 4: Temporal signature of TLS during free-evolution. Echo phase-memory decay 
(dots) measured at e = 256MHz (red), 320MHz (green) and 384MHz (blue), and corresponding 
simulation (solid lines) with assumed flux noise PSD function, S'(f) = A/f 09 + SiC(f;F\,Wi), 
where A = (27r0.9 x 10 6 ) 2 (rad-s -1 ) 2 , S± = 9.0 x 10 8 rad-s _1 , F\ = 1.25MHz, W\ = 0.35MHz. 
The temporal signature of the Lorentzian at 1 MHz is essentially independent of the e bias. Inset 
is the echo pulse sequence. Here, r is the total free evolution time in the sequence. Note that the 
20 MHz bump has little influence at these bias points. 
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Supplementary Figure SI. System Hamiltonian for different reference frames, (a) Lab frame. Static = + £a x , 

Hrf = A r [ cos( 27 T!/ r ff + (f>)a x , and 9 indicates the tilting of quantization axis when transformed into the qubit frame, (b) Qubit 
frame. The red arrow (z/ q ) is the free-evolution quantizing field, and 2ni' r ft indicates the time-dependent transformation into 
the rotating frame, (c) Rotating frame ((f) = 0). The pink arrow (v' R ) indicates the effective driven-evolution quantizing (Rabi) 
field in the rotating frame in the presence of a finite detuning Aia (d) Rotating frame (same as (c) but Au = 0). The yellow 
arrow (v R ) indicates the resonantly driven-evolution quantizing field, or the pseudospin. 
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(a) System dynamics in T, and T t experiments 
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(b) System dynamics in FID and Rabi experiments 
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Supplementary Figure S2. System dynamics in analogous experiments between free and driven evolution, (a) 

7i versus T± p . The figures represent the typical case when the initial qubit state is parallel with the quantizing field (/3 = 0). 
Depolarization shown is based on the condition, i/ R -C fc B T/h<C rz q , so that, in the steady state, (dv) ~ ~ 1 for the Ti process, 
and ( ax) ~ 0 for the T\ p process, (b) FID versus Rabi. The dashed line indicates the precession, while depolarization happens 
simultaneously. The process shown ignores the longitudinal relaxation. 
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Supplementary Table SI. Comparison of decohering properties of a flux qubit during free evolution and driven evolution 
(weak and resonant). 
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SUPPLEMENTARY NOTES 


Our test device is a superconducting flux qubit, an aluminium loop interrupted by four A1-A1CR-A1 Josephson 
junctions. In the presence of external harmonic control, the Hamiltonian of the superconducting flux qubit in the 
laboratory frame (Supplementary Figure Sla) is 


H = - [Ady + ea x + A rf cos(27rzy. f f + (/>)a x \ , 


(SI) 


where a z and a x are Pauli operators, e = 2I p $b/h (/ p is the superconducting persistent loop current) is called the 
energy flux bias, which gives the energy difference between the two classical circulating-current states (clockwise and 
counterclockwise) and is adjusted by the external magnetic flux $ via <Fb = *F — $0/2 ($0 is the superconducting 
flux quantum). A = 5.4GHz is the tunnel-coupling strength between the two states. The oscillating term represents 
the harmonic drive with amplitude A r f, carrier frequency (or microwave frequency) i/ r f and phase <j>. In this device, 
noise in s, A and A r f are considered the main physical sources of decoherence in most experiments, and have simple 
geometry in the lab frame. However, the decohering properties are routinely defined and better understood in the 
frame spanned by the qubit eigenbasis. 

By making the transformation, 


a x , = cos 9 & x — sin 9 a z , 

°y = a y > 

a z > = cos 9 c t z + sin 9 a x , (S2) 

where 9 = arctan(e/A) is the rotation angle of the quantization axis from the lab frame, we can write the Hamiltonian 
in Eq. SI in the qubit (eigenbasis) frame (Supplementary Figure Sib): 

H' = — \v q dy + A r f cos 9 cos(27rzy.f t + <^>)dv + A r { sin 9 cos(2ni/ r {t + <(>)dv] (S3) 

where u q = \/e 2 + A 2 is the qubit’s level-splitting. 

The free-evolution dynamics can be visualized on the Bloch sphere in this qubit frame (Supplementary Figure 
S2). Ignoring the external control, the only term left is the one along the quantization axis. The qubit, like a spin, 
undergoes right-handed (counterclockwise as viewed from the north pole) precession around the ./-axis, at a cycling 
rate equivalent to the qubit frequency. The qubit’s Bloch vector and the z'-axis form an angle /? depending on the 
initial state. For example, if the qubit starts from the ground (| —z ; ), south pole) or excited state (| +z'), north pole), 
/3 = 180° or 0°. Hence the qubit state remains unchanged. When the qubit is initially prepared at a 50%-superposition 
state, e.g., |±a/) or |±z/ / ), /? = 90° and the qubit will precess around the equatorial plane. I 11 general, for an arbitrary 
value of /3, the qubit leaves a cone-like trace which is axially symmetric around z’. 

The free-evolution decohering processes are defined with respect to the quantization axis. Note that, decoherence 
in superconducting qubits should be understood in a sense of time-ordered ensemble average. 

T\ is the longitudinal depolarization (longitudinal relaxation time), or often just “relaxation time” for short. From 
perturbation theory, such longitudinal relaxation is the consequence of environmental modes at the qubit frequency 
transversely coupled to the qubit. The inversion recovery experiment [26], which initializes the qubit at the excited 
state and records that how it temporally relaxes along z' (Supplementary Figure S2a), is one of the standard methods 
for measuring Xj. 

T 2 is the transverse depolarization (transverse relaxation time), or just “coherence/dephasing time”. In the weak 
coupling regime where decoherence of the qubit during the correlation time of the environment is negligible, the 
Bloch-Redfield theory [37-39] gives a simple description of the exponential dephasing rate: I/T 2 = l/2Ti + 1 /T v , 
where T v is the pure dephasing time. Unfortunately, in many instances such as l//-type spectrum, the Bloch-Redfield 
approach does not apply because the noise correlation times become large at low frequency. However, the decay 
functions due to Tj processes and due to T ^ processes still factorize [40]. Pure dephasing is a consequence of ensemble 
averaging over a fluctuating qubit frequency [41]. 

Depending on the experiment, there are various types of “T^. The simplest one is the free induction decay (FID) 
or Ramsey experiment [26], which measures the associated coherence time Tfid with bare precession (one-way phase 
accumulation). The measured pure dephasing is the result of ensemble averaging over a inhomogeneous parameters [14] 
(here, temporal inhomogeneity due to low-frequency noise for superconducting qubits). If a refocusing pulse is added 
in the middle of the free-evolution period as in the spin-echo experiment [36], dephasing due to low-frequency noise 
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can be reduced. In superconducting qubits, the noise PSD is l//-type, and we typically obtain a longer echo phase 
decay time TEcho > Pfid- More advanced sequences such as CPMG [12] (a generalized version of the spin echo 
sequence) may further improve coherence and, in addition, be utilized to extract the free-evolution PSD through its 
noise-shaping/filtering properties. 

Driven-evolution dynamics are conveniently described in a reference frame which rotates around z' at the drive 
frequency zy.f by a second transformation, 

ax = cos(27Tiy.ff)(5v + sin(27rzy.ff)dy , 
ay = cos(27ri/ r ff)<5y — sm(2mv T {t)a x i , 

a z = d y . (S4) 

The Hamiltonian in Eq. S3 is then written as (Supplementary Figure Sic) 

H = — [A v&z + ~H r f cos cos (f> ax + sin ay) 

+ -A r f cos #( cos(—2zr • 2zz r f t — <j>)&x + sin(—2zr • 2z/ r ft — < f>)ay ) 

+ Atfsindcos(2nv T {t + 4>)&z\ , (S5) 

where Azz = u q — zz rf is the frequency detuning between the qubit and the driving field. In the weak driving limit 
(A r t ■Ci'rf), the last two lines in Eq. S5 can be omitted, since these rapid oscillations average to zero on any appreciable 
time scale of the qubit dynamics in the rotating frame (rotating wave approximation). The Hamiltonian within this 
approximation is given by 


H = (h/2) [Au&z + v R (cos<j)ax + shindy)] , (S6) 

where u R = |H r f cos6 is the Rabi frequency under resonant driving (i/ T { = zz q ). In the presence of finite detuning, 
Eq. S6 describes the effective driving field, whose effective Rabi frequency, z/ = u R \J 1 + (Av/iy R ) 2 ss v R + Av 2 /2v R 
is linear with z/ R and quadratic with Aza 

For resonant driving along the X-axis (Au = 0, <f> = 0), the Hamiltonian in Eq. S6 describes a pseudospin quantized 
along X with a level-splitting equivalent to the Rabi frequency v R (Supplementary Figure Sid). Therefore, the driven- 
evolution dynamics represented in the rotating frame are analogous to the free-evolution ones in the qubit frame. That 
is, the pseudospin precesses around the X axis at the rate of z^ R , and its subtended angle depends on the initial state. 
For example, if the qubit is prepared at either |+X) or |—X) state, its Bloch vector will remain in that orientation, 
as does the corresponding free-evolution scenario when the qubit is initialized in its excited (l+^O) or ground (| —z')) 
state. When the initial state is an “equatorial state” (here, lying on the Y — Z plane), the pseudospin will precess 
around the Y—Z plane just as their free-evolution counterparts precess around the x'— y' plane (equatorial states for 
the free-evolution case). 

Following this analogy, we define decohering processes for the pseudospin with respect to the new quantization axis 
X. T\ p is the rotating-frame longitudinal depolarization/relaxation time, i.e. relaxation in the rotating frame, which 
corresponds to depolarization along the X-axis. The environmental noise responsible for this decay is the transversely 
coupled (_L X, comprising fluctuations z^ R of the Rabi-frequency). To measure T\ p , we use the T\ p experiment described 
in this manuscript, which prepares the qubit at the |±X) state, and then records how it depolarizes during continuous 
driving. 

T 2p is the rotating-frame transverse depolarization/relaxation time. Similar to the free-evolution case, this transverse 
decoherence can also be factorized into a Tip-induced part and a pure dephasing part T vp , all with respect to the new 
quantization axis. The rotating-frame pure dephasing is a consequence of fluctuation in the pseudospin frequency. 
It has first-order sensitivity to noise in the Rabi frequency z^ R (due to driving field amplitude fluctuations), which is 
found to be the dominating source of Rabi decay in this device [18]. 

T 2p , like its free-evolution counterpart, also varies in different experiments. The simplest one is the original Rabi 
experiment [26], which measures the decaying oscillatory signal (characteristic time Tn a bi) by resonantly driving the 
qubit from its ground state (the Rabi experiment in a general sense means the driven evolution starting from an 
arbitrary state. Hence, even the T\ p experiment is just a special case of Rabi. However, in this article, we denote 
Rabi as its most original version only, in which the qubit starts from its ground state | —z') (=| —Z))). It is the 
driven-evolution analog of the FID experiment, since both characterize the transverse decay during “bare” equatorial 
precession. The pure dephasing part is sensitive to u R (linear coupling, first order) and Az^ (quadratic coupling, second 
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order) inhomogeneities, both of which modify the pseudospin frequency, i.e. the effective Rabi frequency 1 / . If the 
second half of the Rabi pulse is 180°-phase-shifted as in the so-called rotary-echo experiment [18], the Rabi field in 
the rotating frame points to —A" instead of A, and the precession reverses direction. Therefore, slow fluctuations in 
the pseudospin frequency can be refocused. Rotary-echo is the driven-evolution analog of spin-echo, and the decay 
(characteristic time Tkotary) can be analyzed in a similar way. Furthermore, the analogy to the CPMG sequence is a 
multi-pulse generalization of the rotary-echo sequence. 

The connection between the driven-evolution times (Ti p , ), the free-evolution times (Tf, T 2 ), and the power 
spectral density (PSD) of the physical noise terms (noise in A and e) is formally made using the generalized Bloch 
equations (GBE) [23, 24] and the transformations between Eqs. S1-S6. 

The qubit-bath interaction Hamiltonian is defined as 

1-ii = ^ JQa g z + Qe g x ] , (ST) 

where Qa and Q e represent the bath variables which couple to the effective tunnel-coupling and flux terms, respec¬ 
tively. 

When transformed into the qubit frame, the Hamiltonian becomes 

H'l = 7 ^ [(Qa cos 9 + Q e sin #)dy + (-Qa sin 9 + Q s cos 9)<j X '] . (S 8 ) 

Therefore, the PSD in the qubit frame is a (co)sinusoidally weighted sum of the PSD in the lab frame, namely, 

S x {f) = sm 2 9S A (f) + cos 2 9S^f) , 

SAf) = cos 2 e S A (f) + sm 2 9 S £ (/) . (S9) 

S\(f) stands for the PSD of A-noise. It is defined as the reduced power spectral density, i.e., S\(f) = (27r) 2 dti|(A(0)A(t)+ 

A(t)A(0)) exp(— i2nft) and |(A(0)A(t) + A(t)A(0)) is the symmetrized autocorrelation function. In the classical regime 
(hf < k B T , where T is the effective temperature), it can be replaced by the classical autocorrelation function, 
(A(O)A(t)). A G {x',y', z'... or A(=Q A ),e(=Q e )... } indicates either the fluctuating part of system Hamiltonian or 
the axis along which those fluctuations occur. 

I 11 the rotating frame, the stochastic Hamiltonian becomes 

'Hi -~ ^ [(Qa cos 6 + Q e sin 9)a z + • • • 

+ (—QAsin# + Q e cos 9) cos(2Tns T {t)&x — 

— ( — Qa sin0 + Q e cos9) sin(27rzz r ft)<jy] . (S10) 

This Hamiltonian, together with the deterministic one in Eq. S 6 , governs the system’s motion for a single realization. 
Comparing Eq. S 8 and Eq. S10, ^/-noise in the qubit frame is invariably transferred into Z-noise, while aCnoise is 
separated into A-noise and F-noise with a frequency mixing of iy.f, indicating that the rotating-frame noise Hx,y{f ) 
is actually modulated (frequency-shifted) qubit-frame noise, derived from S x >(\f ± v v {|). 

In the weak coupling regime, ensemble averaging over all the heat bath realizations gives the non-Markovian 
equations of motion, from which the GBE are derived. With resonant driving and rotating wave approximation, we 
have the simplified GBE: 


d(dy) 

dt 

d(dy) 

dt 

d (&z) 
dt 


-Tx{g x ) + v x , 

-Py((Ty) + lS K (a Z ) + V Y , 
—r z(gz) ~ V r (<7y) + Vz ■ 


The depolarization coefficients are determined by the PSDs of the heat bath fluctuations: 


r * = \ [Sx{v ,q + Ci) + S x >(v q ~v n )]+\ . 

r Y = + 

r z = ^ S X ’(yA + — \S x ’(y q + i^) + S x ’(y q — zz R )] , 


(Sll) 


(S12) 



and v\ (A = X, Y, Z) (see [23, 24] for details) defines the steady-state polarization: 


{&x 


SS 


<^r 

<^> ss 


Vx 

r x ’ 

T z v y + v R v z 

+ rvr z 

V Y v z - v R v Y 

vl + r Y r z 


(S13) 


If the system-environment coupling is weak enough such that the driven-evolution depolarization is negligible on the 
time scale of qubit precession in the rotating frame, i.e. v R (A = X,Y,Z), we have (& Y ) SS = (d\z) ss = 0. In 

addition, (&x) ss = 0, if v R <^k B T/h. 

When there is finite detuning (Azz 7 ^ 0), the quantization axis of the pseudospin becomes the effective Rabi field, 
i.e. the pink arrow in Supplementary Figure Sic. I 11 this case, the longitudinal steady-state polarization is expressed 
by 


sin ?7 S x {v q ) tanh() + cos 2 /? S z {v R ) tanh( 
\ (1 + sin 2 ^) S x {u q ) + cos 2 7 / S z {v R ) 


hU-r, 


{cr Z ') = 


(S14) 


where Z' indicates the longitudinal axis of the pseudospin in the presence of small detuning (A v < v R ) 1 and zy = 
arctan(Az^/T R ) is its tilted angle from the equatorial plane. In both the weak driving and weak coupling (driven- 
evolution) limits, and to first order in detuning (or rf) 1 the steady-state X-polarization is 


(dx) ss = COS V (& z .r 


smr)S X f(is q ) 
l Sx{v q) + S z {v R ) ' 


(S15) 


The component (ax') ss is negligible in the weak coupling limit. 

Connecting the results in Eq. S12 with the T\ p and Rabi experiments, we have 

Tip = Tx and = F^ abi = — (Ty + T z ) = — Fi p + T vp , (S16) 

where Y vp = | S x '{v q ) is the rotating-frame pure dephasing rate in the Rabi experiment. Note that, in this model, 
we did not take into account the low-frequency fluctuations in Az^ and v R , which, in fact, dominate the Rabi decay. 
Therefore, we use the “f” symbol to indicate that the rate is free from those inhomogeneities. However, decay due to 
Az^ and v R inhomogeneities can be calculated semiclassically, in the same way as is done in the FID case [41]. 

As previously discussed, the longitudinal relaxation rate Fi p depends on the transversely coupled (rotating-frame) 
noise at the frequency of v R , i.e. Sz(is R ) and Sy{z/ R ). While Sz{v R ) is equivalent to S z '(v R ), Sy{y R ) is derived from 
S x '(v q ± u R ) due to the frequency mixing during frame transformation. This explains all the PSD components that 
contribute to Ti p . On the other hand, the transverse decoherence r^ abi combines a relaxation-induced part 5 Tip 
and pure dephasing T vp in the rotating frame. The presence of the qubit-frequency noise S x >(u q ) in the expression for 
T vp is due to the fact that, when transferred into the rotating frame, S x '(v q ) is mixed down to zero-frequency noise, 
S'x(0), which acts as the longitudinal quasistatic noise leading to dephasing. From above, we see that the decoherence 
formalism of driven evolution in the rotating frame is truly analogous to that of free evolution in the qubit frame. 

At weak driving, we have v q ±i' R sst'q, which further simplifies Eq. S16: 

Fi p = - T], + r„, 

F Rabi = ( S17 ) 


where 


Ti = ^S x {v q ) = i [sin 2 #SA^q) + cos 2 #^(z/q)] , 
= \ Sz'i ” R ) = \ [cos 2 dS A (v R ) + sin 2 0 S^ R )\ . 


(S18) 
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ri is the free-evolution longitudinal relaxation rate, and T u is the rate associated with the Rabi-frequency noise. Their 
relation to Ti p inspires us to extract noise PSD at the Rabi frequency by measuring the lab-frame and rotating-frame 
longitudinal relaxation rates and taking the subtraction. 
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